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Abstract 

We show that quantum optical systems preserving the total number 
of excitations admit a simple classification of possible resonant transi- 
tions (including effective), which can be classified by analizying the free 
Hamiltonian and the corresponding integrals of motion. Quantum sys- 
tems not preserving the total number of excitations do not admit such a 
simple classification, so that an explicit form of the effective Hamiltonian 
is needed to specify the allowed resonances. The structure of the resonant 
transitions essentially depends on the algebraic propereties of interacting 
subsystems. 

1 Introduction 

In a direct analogy with classical mechanics, composed systems in quantum 
optics (which describe interaction between several subsystems) can be divided 
into two classes: 1) Systems which possess a necessary number of integrals of 
motion, so that, the classical counterpart is an integrable system; 2) Systems 
that do not admit a sufficient number of integrals of motions, so that, the 
classical counterpart is a non-integrable system pQ . A quantum system can have 
basically two types of integrals of motion: a) Kinematic integrals, which do not 
depend on the kind of interaction between subsystems, as for instance, the total 
number of atoms; b) Dynamic integrals which arc related to the particular form 
of interaction between subsystems, as, for instance, a number of excitations 
preserved in some transitions between the energy levels of the subsystems. 

Typical for quantum optical systems dipole-like interactions between two 
subsystems (X and Y) can be described with a generic multichannel Hamilto- 
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nian of the following form: 

H = J2 "j* 3 * + E + E 9* (xi + x{) + Y 



(1) 



where the two first terms represent the free Hamiltonians of the subsystems, 
so that the frequencies u)j,£lk > 0, and the last term describes the interaction 
between them. The operators Xq ± , Yq ± , [Xq ±,Yq ±] = 0, are usually elements 
of some deformed algebra [2], and in particular, satisfy the ladder commutation 
relations 



x ,x± 



±X> ± , 



±Y£. 



(2) 



In the interaction Hamiltonian there are two kinds of terms: of the form 
X 3 _Y+ and Xl_Y_. It is easy to observe that in the rotating frame, that is, 
applying the following unitary transformation 



U = cxp it 



E^ x o 



k 



n k Y* 



to the Hamiltonian |JI]|. the counterrotating terms ~ X-^Y^ oscillate in time with 
a frequency u)j + fife and the rotating terms ~ X J _Y+ oscillate with a frequency 
ujj — Qk- It is clear that under the condition Uj w fife the rotating term XiYf 
in Q is approximately time independent (and thus, can generate transitions 
with a probability of one between the energy levels of the system), meanwhile 
the counterrotating term, XlSY^ always oscillates rapidly in the rotating frame, 
and its temporal average is zero. 

By neglecting the counterrotating terms in the Hamiltonian (|TJ, which is 
commonly called the Rotating Wave Approximation (RWA), we arrive at the 
Hamiltonian 



H = Y, "jX J o + E ^ Y o + E ft* { X - Y + + X + Y - 

j k j,k 



(3) 



which admits several dynamic integrals of motion iV p , generally not allowed in 
|QJ. This implies that the whole representation space of the system is divided 
into finite dimensional invariant subspaces, and the mathematical treatment is 
essentially simplified. 

It is worth noting, that the semiclassical models, when some of the sub- 
systems are described by c-numbers instead of operators, are treated using es- 
sentially the same type of Hamiltonians as l|T]l. The semiclassical transition 
in Q can be done by going to the rotating frame of the semiclassical system 
and then, just substituting the transition operators by some complex numbers. 
For instance, in the case of single channel interaction, j = 1, k = 1, when the 
system Y acquires classical features, the Hamiltonian in the rotating frame 
corresponding to the system Y takes the form 



H = ujX + g (A_ + X+) (F_e 



iClt 



■Y, 



-i(lt\ 
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so that the corresponding semiclassical Hamiltonian is obtained by substituting 
Y+ -> 9,Y- -> g*, giving 

H sc = loXo + g (X_ + X+) (g*e mt + ge- mt ) . (4) 

Such Hamiltonians usually appear when a quantum oscillator and/or a collection 
of atoms is pumped by an external force [3], @], 

In the RWA-like systems, described by Hamiltonians of the form (0), the 
resonance conditions 

Uj « Q k , (5) 

means that the term X 3 _YK which explicitly appears in the Hamiltonian, does 
not depend on time in an appropriate rotating frame. Nevertheless, such ex- 
plicit resonances are not the only kind of resonances, which can be found in the 
Hamiltonian J3J . Usually, the composed systems admit several types of implicit 
resonances related to effective transitions, which do not appear in the original 
Hamiltonian, between their energy levels. Such effective interactions play im- 
portant roles in many physical applications and can be revealed by adiabatic 
elimination of slow transitions [Hj- The implicit (effective) resonances are char- 
acterized by their position, • mjCOj ~ ^2 k nk^lk, and strength, i.e. in what 
order of perturbation expansion they appear for the first time. Although, a 
generic system can possess a large number of different types of effective transi- 
tions, all the possible resonance conditions can be classified only by analyzing 
the free Hamiltonian and the integrals of motion. Because all the invariant 
subspaces are finite-dimensional, there are always a finite number of different 
resonances. We will refer to these kinds of resonances as kinematic resonances, 
which include both explicit and implicit resonances. In Sec. II we show with the 
example of atom-field interactions, that it is possible to classify all the kinematic 
resonances in a straightforward way. 

The situation is quite different in quantum systems with a lack of integrals 
of motion, corresponding to classically non-integrable dynamic systems. In such 
systems an infinite number of different resonances arise and a priori it is im- 
possible to determine their position and strength, which essentially depend not 
only on the type of interaction but also on the algebraic properties of each inter- 
acting subsystem 0, jS], In Sec. Ill we will discuss such dynamic resonances 
analyzing different models of interaction of quantum and classical fields with 
atomic systems. 

2 Kinematic resonances 
2.1 A simple model 

An important example of kinematic resonances is the Dicke Model ^Uji this 
model describes the interaction of a collection of A identical two-level atoms 
with a single mode of a quantized field under the Rotating Wave Approximation. 
The Hamiltonian that governs this system is given by 

H = ujo) a + ujqS z + g(aS+ + a'S-), (6) 
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A 

where to is the held frequency, loq is the atomic frequency, and S z = ^2 Szi , S± — 

i=l 

A 

Y] s±i are the collective atomic operators, they represent the atomic inversion, 

i=l 

and the transition between the atomic energy levels respectively, their commu- 
tation relations are given by the su(2) algebra, 

[S Z ,S±] = S ± , {S+,S-]=2S Z , (7) 

and , a are the creation-annihilation held operators, obeying the bosonic com- 
mutation relations, [a,a^] = 1. 

The Hamiltonian JBJ admits two integrals of motion, a kinematic integral 
of motion, given by the total number of atoms (that is constant for a closed 
system), A = S 11 + S 22 , where S vl is the atomic population operator for the 
i-th atomic energy level, and the dynamic integral of motion, corresponding to 
the total number of excitations, has the form 

N = S z + a) a, [N, H] = 0. (8) 

This system is the simplest non trivial example of one-channel quantum 
transitions: an absorption of one photon is accompanied by an excitation of one 
atomic transition. Using the integral of motion (JHJ , the Hamiltonian @ can be 
rewritten as follows, 

H = uN + H int , 
where the interaction Hamiltonian is 

H int = AS z + g(aS + + JS-), (9) 

and A = luq — w, is the detuning between the held and the atomic transition fre- 
quencies. If the only possible resonant condition, to = ujq, is held, the interaction 
Hamiltonian © is reduced to, 

H int =g(aS + + d*S-), (10) 

which implies that the atomic transition probability (as a function of time) 
oscillates between zero and one. 

On the other hand, in the far-off resonant (dispersive) limit, |A| 3> gV^, 
the interaction Hamiltonian 10 is diagonal, 

H mt ~AS Z + ^ [(2a f a + 1)5, - S 2 ] , (11) 

so that only the phase of the system evolves due to the appearance of the so- 
called dynamic Stark shift terms ^21 m the effective Hamiltonian Qllfl. 
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2.2 Generic atom-field interactions 



Let us consider the interaction of a system of A identical iV-level atoms of an 
arbitrary configuration with a single mode of a quantized field of frequency 10. 
The Hamiltonian describing this system has the form 

H = H d + H int , (12) 

where 

N 

H d = uJa + E i s " > E 3 < E i+i > ( 13 ) 

3=1 

is the free Hamiltonian, and S n are the collective atomic population operators 
corresponding to the j-th atomic level of energy Ej , and Hint is the interaction 
atom-field Hamiltonian, whose explicit form depends on the atomic configura- 
tion, in the dipole approximation, i.e. only one-photon transitions are allowed. 
In this Section we suppose that the Rotating Wave Approximation is imposed, 
so that the total number of excitations in the atom-field system is preserved, 
and thus, the whole representation space of this quantum system is divided into 
finite-dimensional invariant subspaces. 

The Hamiltonian l|12|) admits two integrals of motion: a kinematic integral, 
given by the total number of atoms: 

N 

A = Y / S l \ (14) 

i=l 

and a dynamic integral, corresponding to the total number of excitations in the 
system: 

N 

N = a^a + J^^S 11 , (15) 

i=i 

where the parameters \ij depend on the atomic configuration. 

Let us note that a generic interaction term can be written as follows 

at-1 k 

f(S u ,aU) J] (Sf ) V", l. [r /. .;, ! V. (16) 

3=1 

where S % ± (SV — , S 1 ^ = j > i) are the atomic transition operators 
satisfying the u(N) commutation relations, [S v , S km ] — 8i m S k i — 5kjS lm , and 
negative exponents correspond to the Hcrmitian conjugated operators. The 
operational coefficients f(S li , a) a) depend on diagonal atomic operators and the 
photon number operator. In what follows we will omit the coefficient f(S 11 , a^a), 
since it leads only to some phase shifts and does not change the distribution of 
excitations in the system's energy levels. It is worth noting that (|16fl is not a 
unique way to represent a generic interaction term. 
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Since the Rotating Wave Approximation is imposed, every interaction term 
should preserve the total number of excitations, so the condition 



N. 



JV-l 

n 

3=1 



5: 



-.few 



is held, and thus, the numbers kj satisfy the following restriction 

JV-l 

k N = ^2 k j (a*JV - H) • 
Thus, any admissible interaction can be described as 



JV-l 

n(«? 



JY 



N — l 



(17) 



(18) 



The interaction (|18[) becomes resonant when the atomic transition energies and 
the field frequency satisfy the following condition 



JV-l 



^ kj (E N -Ej-w (fiN - Mi)) = 0- 



(19) 



i=i 



The important point here is that interaction terms l |18| l can be explicitly pre- 
sented in (|12|) . or can describe effective interactions, and thus should be obtained 
from the original Hamiltonian by adiabatic elimination of some far-off resonant 
transitions. So that, the condition (|T§1) describes explicit resonances if the cor- 
responding interaction term is present in the original Hamiltonian or implicit 
resonances, if such interaction is effective. It is clear that the total number of 
both explicit and implicit resonances is finite, which is a consequence of the 
restriction (|17|) imposed by the Rotating Wave Approximation. 

The number of possible resonances depends on the number of atomic levels, 
N, and the total number of atoms, A. Note that the resonance condition (|19fl is 
associated with a vector k = (k±, ...fczv_i), in order to not repeat the resonances, 
we consider only vectors k with coprime components (that is fci, ...fc/v-i do not 
have a common factor). Then, different vectors k satisfying the following 
condition 



AT-l 



JV-l 



i •■••) l^*li=l,...,JV-l ^ < ^ 



(20) 



i=l,...,JV-l 

define different resonances. 

On the other hand, taking into account the dynamic integral of motion (|15(l . 
one can rewrite the free Hamiltonian Ijl3|l as follows, 



JV 



H d = ujN + {Ej - Hjw) S". 



(21) 
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Summing the kinematic integral —j3 J2i=i to the above Hamiltonian, 

where the constants (3 = (3(n), n = (m, .., Jijy)) are chosen such that, 

N N 

= n 'J ( E i ~ ^3^) , ^ Tij = 1, nj-€2, 

we obtain 

JV 

2? d = wJV + ^ (2^- - /xj-w - /3(n)) . (22) 

3=1 

It is easy to see that the condition Ej — [ijU) — f3(n) = 0, for any fixed j and 
some values of (3(n), enumerates all the possible resonances (|19|l . Let us take 
for instance j = 1, if Ei — /Ltiw — /3(n) = 0, then 

2V-1 

^ rij (E N - Ej —uj (fiff - fij)) + (ni - 1) (E N -Et-w (pn - Mi)) = 0, 

(23) 

which coincides with (|19|) . when fcj = for j = 2, ...JV — 1 and fci = (ni — 1). 

This means that we can always represent the free Hamiltonian in a way 
that all the possible resonance conditions, corresponding to both explicit and 
implicit resonances, appear as zeros of coefficients of the atomic population 
operators S 33 in the free Hamiltonian, after taking out the integral of motion 
l|15|) corresponding to the total number of excitations. 

The above allows us to classify all the possible kinematic resonances: 

1. Multiphoton resonances: transitions which involve absorption and emis- 
sion of photons: a) simple n-photon transitions, described by the terms 
~ a n S k ^ , with the resonance condition Ej — Ek ~ nuj. The terms with 
n = 1 can be present in the original Hamiltonian and in such case rep- 
resent explicit resonances. Multiphoton transitions appear even in a sin- 
gle atom case; b) Collective atomic transitions, described by the terms 
~ a n SfS^, (in general a product of several, up to A atomic transi- 
tion operators can appear), with the corresponding resonance condition 
(Ei — Ek) — (Ei — Ej) rs ncu. It is clear that such resonances can appear 
only in multi-atom systems. Note, that if I = i (or k = j) then we obtain 
the same resonance condition as in a), describing an effective process of 
absorption of n photons with atomic transition from fc-th to j-th energy 
levels. Nevertheless, the corresponding term in the effective Hamiltonian 
would be multiplied by the atomic population operator S u , which means 
that more than one atom is need to realize such a process. 

2. Virtual photon resonances: atomic transitions between independent chan- 
nels caused by the quantum field fluctuations, and thus existing even when 
the field is in the vacuum state. Such resonances appear only when the 
system has more than one atom and are described by terms in the form of 
products of atomic transitions operators. For instance, the simplest term 
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of this kind (typically appearing in the first order perturbation expansion) 
is S 1 !, k ^ i, j 7^ I, describes atomic transitions k — > j, I — > i, and the 
corresponding resonance condition is 

{Ej - E k ) - (Ei - Ei) « 0. (24) 

Obviously, such transitions can be realized in a system which consists of 
at least two atoms. Note that the term ~ S+S 1 !*, which represents the 
atomic transition k — * i, I — > j, satisfies the same resonant condition. More 
involved interactions, like~ S+S^S™ n , can appear in the highest orders 
of the perturbation theory. The strength of virtual photon transitions does 
not depend on the field intensity in the leading order of the perturbation 
expansion. 

3. Photon assisted transitions: atomic transitions when every photon emis- 
sion is accompanied by a simultaneous photon absorption. These transi- 
tions appear only when the atomic configuration contains coherent chan- 
nels, similar to lambda-like configurations. The strength of such interac- 
tions depends on the number of photons in the field and the populations 
of some atomic energy levels. The simplest terms describing the photon 
assisted transition (which appear in the lowest order of the perturbation 
theory) have the form: f{a^a, S" )S+ , corresponding to the resonance con- 
dition Ei w Ek, where the transition k <-» I is not present in the original 
Hamiltonian and f(ata, S^) is a linear polynomial of the photon number 
operator and the atomic population operators. 

In the above classification we do not consider interactions corresponding 
to powers of terms describing some interactions. For instance, the n-photon 
transition corresponding to the term (aS^) n , that describes an absorption of n 
photons by at least n atoms in the transition j — > k, we include in "one photon 
transitions" . 

In the case of interaction of an atomic system with classical field the Rotating 
Wave Approximation implies that the interaction Hamiltonian can always be 
reduced to a time-independent form. In this case, apart from explicit resonances, 
several types of implicit (effective) resonant transitions can take place. It is clear 
that no interactions similar to virtual photon resonances can arise. Nevertheless, 
transitions similar to multiphoton resonances ~ ,k — j > 2 and photon 
assisted resonances ~ S J + , where the transition k <-> j docs not exist in the 
initial Hamiltonian, actually appear. 

2.3 Four level diamond configuration atoms 

As an example of the kinematic resonances we will study the interaction of A 
four level diamond configuration atoms with a single mode of a quantized field 
under RWA. The Hamiltonian governing the evolution of this system has the 
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form 

4 

H = oja)a + E o S " + 5i ( aS + + a ^ S - 2 ) + -92 (aSf + a^S 13 ) (25) 

3=1 

+ 53 {aSf + a^ 24 ) + 9i (aSf + cSS?) . 

This Hamiltonian describes four one-photon atomic transitions, (1 <-» 2, 1 <-> 
3), (2 <-> 4, 3 <-* 4), gathered in two pairs of coherent quantum channels. 
The corresponding resonance conditions (explicit resonances) are E% — Ex ~ 
u, Es — Ex ~ u, E4 — Ei f» w, and E4 — E3 ~ uj. To find all the other 
(implicit) resonances we will follow the method outlined in previous Subection. 
The dynamic integral of motion for this system is 

N = Ja + S u -S 11 , 

so that the coefficients /ij in Q15[l. are fj,i = —1, fi 2 — fJ-3 = 0, ^4 = 1. Substi- 
tuting these /Ltj into l|19f) we obtain all the possible resonance conditions: 

fci (£ 4 - J5i - 2cj) + fc 2 (£J 4 - E 2 - ui) + k 3 (E 4 -E 3 -cj)k0, (26) 

kx,k 2 ,k 3 € Z, and the corresponding (effective) interactions are 

(5| 4 ) fel (5| 4 ) fc2 (5+ 4 )" 3 a *3+fc2+2fc!_ 

The vector fc = (ki,k 2 ,k 3 ) associated with possible resonances should satisfy 
the condition (|20|) 



max • 



+ ^1^,1^1 ^ < A i,i = 1,2,3. 



For a single atom we have the vectors (1,0,0), (0,1,0), (0,0,1), (1,-1,0), 
(1,0,-1), (0,1,-1), which correspond to the following resonance conditions 
(interactions): £4 - Ex « 2w (a 2 ^ 14 ), £4 - £2 ~ w (aS 24 ), E^ - E 3 ^ u 
(aSf) ,E 2 ~Ex~lo (aSl 2 ) ,E 3 -Ex^uj (a^ 3 ) , £3 ~ £2 (Sf) , respectively. 
For two atoms there there are 15 other vector apart from 6 shown above. 

For a better understanding of the nature of the implicit resonances, we find 
the first order effective Hamiltonian considering that all the interactions ap- 
pearing in the initial Hamiltonian (|25|l are far from resonance. Following the 
method outlined in the Appendix, we apply to the Hamiltonian (|25l) the follow- 
ing sequence of unitary transformations 

H eff = U A U z U 2 UxHU\ululu\, 

where Ui = exp (siAi) , 

Ax = aS^-a^S 12 , A 2 = aS™ - a^S 13 , 
A 3 = aSf-a^S 2 *, A4 = aSf - a* S 3 _ A , 
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Figure 1: Interactions diagram for a single four level atom of diamond config- 
uration. Solid line: explicit resonances. Dashed line: Implicit resonances: two 
photon transition, photon assisted transition. 



and the small parameters Ej <C 1 are given by 

_ 9i _ 92 

£ 1 — T, F 7 £ 2 — 



E2 — Ei — uj E3 — E\ — oj 

* - .93 o 94, 

£ 3 — 7; F ' e 4 



i?4 — i?2 — ^ -^4 — -^3 — ^ 

we obtain, in the first order on small parameters, the following effective Hamil- 
tonian 

4 

H = Loa^a + ^EjS^ + <P{a^a 1 S J: >)+g2e 1 [Sf (S 11 +a f a + l) + /i.e.] 
i=i 

+ 54 e 3 [S 23 - «M + ^- C -] " (03£i + 54£ 2 ) (a 2 5| 4 + fc.c.) 
+ 53 £i (SfSi 2 + /i.e.) + .g 3 £ 2 (SfSi 3 + /i.e.) 

+ 54 £i (S^S 34 + h.c.) + g A e 2 (SfS 34 + h.c.) . (27) 

where 

$(ata, fl*) = <?i£i (S 22 (S 11 + 1) + ot a (S 22 - S 11 )) 
+92S2 (S 33 (S 11 + 1) + o+a (S 33 - S 11 )) 
+9se 3 (S 44 (S 22 + l)+ala(S^-S 22 )) 
+g 4 e 4 (S M (S 33 + 1) + ota (S 44 - S 33 )) 

is the dynamic Stark shift |11|. 

Let us classify the effective resonances present in the above (first order)Hamiltonian: 
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1. Multiphoton resonances: The only two resonances of this type are one- 
photon and two-photon resonances. There are four one-photon transitions 
corresponding to explicit resonances, which do not appear in the effective 
Hamiltonian J2ZJ- The term ~ a 2 S] 4 describes two-photon transitions, 
with corresponding resonance condition E± — E\ ss 2u>. 

2. Virtual photon resonances: the last four terms in (|27|l . ~ S^S 12 , S+S}?, 
S^S 34 , SfS 34 -, with the resonance conditions E 4 - E 2 « E 2 - E x , E A - 
E 2 ~ £3 — Ei , E 2 — Ei w £4 — _E 3 , and E3 — E x & E^ — E 3 respectively. 

3. Photon assisted resonances: The only resonance of this type is of the form 
~ f{S 4 \S 1 \a^a)S 23 , is produced in the first order of the perturbation 
theory generating transitions between the middle atomic levels, and the 
corresponding resonance condition is E3 ps E 2 . 

3 Dynamic Resonances 

In this Section we study quantum systems corresponding to classically non- 
integrable systems, mainly focusing on the models describing interaction of 
atoms with quantized and classical fields without applying the Rotation Wave 
Approximation 0, [HI (some different systems not preserving the number of 
excitations are discussed in |SJ). 

It is convenient to start with a general analysis of such systems. The main 
idea consists in removing the counter-rotating terms from the multi-channel 
Hamiltonian using a method of adiabatic elimination. In the limit of weak 
interaction, gjk <C uj n ,Q m one can apply, for instance, the Lie-transformation 
method outlined in the Appendix. Although, such analysis can be performed 
for a general system, we will focus on the simplest case of a one channel Hamil- 
tonian and show, that even in this case the absence of the integral of motion 
corresponding to the total number of excitations leads to the appearance of a 
series of (dynamic) resonances, which can be classified according the interaction 
type. 

Consider a single-channel Hamiltonian |{T{| 

H = ujX + flY + g(X + Y_ + X_Y+) + g {X+Y+ + , (28) 

where Xq (Yq) is the free Hamiltonian of the subsystem X (Y) and X + (Y + ). 
X- (Y-) are the up and down operators (see Appendix), respectively, which 
describe transitions between the energy levels in the subsystem X (Y), they 
hold the su c i{2) 2 commutation relations From now on, we do not impose 
any commutation relation between the transition operators, which are generally 
some functions of the diagonal operators and integrals of motion [iVi,Xo,±] = 



[Jv 2 ,y ,±] = 0, 



Pi(X )=V Xo cl>i (X ,Ni), 
P 2 (io)=Vr o 2 (Y ,N 2 ). 
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where Vxo&i PQ), -Ni), Vy ^ ftj ^2) are given by the structural functions 
(appendix (|49|l ). and, in general, (we are going to omit the dependence of the 
integral of motion) 

v,0(*) = 0O?)-0(« + i). 

The consequences of the existence of the counterrotating term ~ XY in the 
Hamiltonian (|28l) are: a) the dimension of the whole representation space is 
the product of the corresponding dimensions of the subsystems: dim(A ® Y) = 
dim A" • dimY", b) there are some additional resonances apart from u) = SI, 
that in the case of a single-channel Hamiltonian do not exist, c) The type 
of resonances depends on the structure of algebras describing the X and Y 
systems. 

The counterrotating term rapidly oscillates (with frequency lo+SI), and thus, 
can be eliminated by applying the Lie-type transformation 

U = exp [e (X+Y+ - , e = g] (u + SI) < 1, (29) 

to the Hamiltonian l|28|) and using the standard perturbative expansion, sec 
Appendix . From now on we suppose that u > SI. 

The elimination of the above term leads to the appearance of new elements 
in the transformed Hamiltonian. All these new terms can be divided into three 
groups: the first group contains the non-resonant terms, of the form X+YV 1 + 
h.c., that can always be eliminated under the condition lu,SI 3> g by applying 
some suitable transformations; the second group consists of resonant terms, 
that cannot be removed if certain relations between u and SI hold, since the 
interaction becomes resonant ( the transformation which eliminates a given 
term from the Hamiltonian becomes singular). This group contains terms of 
the form X+Y™ + h.c., which describe transitions between energy levels of the 
whole system. The third group includes the diagonal terms (functions only of 
X ,Y ), that can never be removed. Our strategy consists of keeping in the 
Hamiltonian only diagonal terms and resonant terms, and, we conserve only the 
leading order coefficients in these terms. 

All the counterrotating terms, like ~ A"Y7™, can be eliminated by applying 
the Lie-type transformation 

U = exp [*(ei..e k )ri (X"Y_J™ - X^Y m )] , (30) 

where n sa gj (jiuj + mSl) <C 1, and 7r(ei..£fc) is the product of powers of some 
small parameters Sj « 1, j = 1, k. In fact %{e\..£^) is proportional to the 
coefficient of the term ~ X"Y" 1 + X™Y™. Applying a sequence of appropriate 
transformations H30J1 we obtain the effective Hamiltonian [J] 

H eff « coX + SlYo + 5 e$ (X , Y ,s) + 



00 00 (_X\l+k-lJ+2(k-l) 

+ ggE (fc _ !)■ {l + k i)j [* k + Y- l+k ^ -Y... K,- --) ! h.r. . CJI, 



k=i 1=0 
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where S = g/2Q < 1. 

The term $ (Xq,Yq,e) represents the dynamic Stark shift (or Bloch-Siegert 
shift) [33] and can be expanded on powers of the small parameter e as follows 

$ (X Q , Y , e) = 4>! (X Q ) cf> 2 (Y ) - 0x (X + 1) 2 (F O + 1) + 0(e 2 ). (32) 

The terms of the form XlfY_ + h.c. describe all the admissible resonant in- 
teractions. In particular the term with k = 1, I = represents the principal 
resonance. The coupling constants 8ki(X , Yq, e) depend on the algebraic struc- 
ture (|49f) of the operators describing both subsystems (more precisely, on the 
degree of the structural polynomials <fii (X ,Ni) and 4>2 (Yo,^)), and in par- 
ticular, they can become zero for some k, I The coefficient in the principal 
term has the form 

9 10 (X ,Y ) = 1, 

and 

i+[fc/2] 

9h(Xo, Y ) = ]T C L[k/2] ( p i ( X o + k)) l+[k/2] - J (-Pi (X )y R k (X Q , Y -2j), 

3=0 

(33) 

R k (X , Y ) = £ CjU, (Pi (X + k))^- 1 (-Pi (X Q )Y V^V 2 (Y - 2j) 

3=0 2 
fe-2(j'+l) 

x n ^ (Y + i) , 

i=l-2j 

for Z > 1, where C 3 k are the binomial coefficients and V°/(z) = 1. Note that the 
product in the last equation is equal to unity if the upper limit is less than the 
lower one. 

In contrast to the case of kinematic resonances discussed in the previous Sec- 
tion, the number of possible resonant interactions appearing in the Hamiltonian 
l|31|) is infinite. These interactions [dynamic resonances) can be classified as fol- 
lows: a) Principal (explicit) resonance, cj = f2, corresponding to the interaction 
term that is explicitly present in the original Hamiltonian ~ X + Y- \ b) Higher- 
order resonances: (21 + k) £1 = koj where I, k G N (I > 1, k > 1), corresponding 
to the effective interactions ~ X^_Y^ l+k , which can be divided into: b.l) odd res- 
onances: when k = 1, 1 = 1,2, 3, that is u> — (21 + 1) f2, b.2) even resonances: 
k = 2, 1 = 1, 3, 5..., with the resonance condition u> = 2mQ, m — 1, 2, b.3) 
fractional resonances: u> = (21 /k + 1) fi, where k and I are coprime numbers. 

It is worth noting that in the vicinity of each resonance (21 + k) O = kuj, only 
the interaction term ~ jffey2/+fc surv i vegp This means that if some resonance 
condition is held, the system is in an approximate invariant subspace and there 
exists an approximate integral of motion N — (21 + k) Xq + kYg. 
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This simple example shows that even in one-channel systems an infinite 
number of resonant interactions may arise, if the total number of excitations 
is not preserved. The appearance of such dynamic resonances in chaotic-like 
systems is expected from the point of view of classical dynamic systems pQ. 
Nevertheless, the quantum nature of interacting subsystems imposes certain 
restrictions on the possibility of surviving of the dynamic resonances. We will 
discuss such restrictions for two examples of the interaction of an atomic system 
with quantum and classical fields. 



3.1 Atom-quantized field interaction 

Let us consider a collection of A identical two-level atoms interacting with a sin- 
gle mode of a quantized field (the Dicke model) without RWA. The Hamiltonian 
that describes this system is 

H = Aa f a + luS z + g(aS+ + a f S-) + .g(a t S' + + aS-), (34) 

and the condition to, A 3> g is held. 
The following identifications 

X = S Z ,X± = S±, Y =n,Y+ = a\Y- = a, 

so that MX ) = (1 + A/2)A/2 - X 2 + X , fc(Y ) = Y , lead to 9 kl {S z ,n) = 
0, k > 3 and the effective Hamiltonian l|31l) takes the form 

H eff sa LoS z +Sla)a + ge{S 2 z + (2a t a + l) S z - {1 + A/2)A/2) (35) 

(=0 ' m=l ^ 

where e = g/ (A + oj) , S = g/2Q < 1. 

Note that the Hamiltonian l|35|l contains only the principal (explicit) reso- 
nance (I = 0), the even and the odd order interactions, but no fractional reso- 
nances. This happens because the structural functions for our subsystems are 
of a first degree polynomial of the photon number operators a^a and a second 
degree polynomial of the atomic population operator S z . 

The integrals of motion for the series of odd and even (exact) resonances are 

N (ki) = (2l + k)S z + ka' i a + 2k5(S + + S-)( a' + a) — 
2k 2 S 2 

-iaTkj [ {21 + k) Sz (2ata + kS ^ + ° ^ ' 

where I = 2m + 1 in N^ 2l \ m = 0,1,2,3,.. and k — 1,2 for even and odd 
resonances correspondingly. 

Let us recall that the higher resonances appear in the effective Hamiltonian 
l|35|) only under the condition oj > A. In the opposite case, u> < A, only the 
principal resonance survives and in the approximation l|31f) the whole effect of 
the counter-rotating terms reduces to the dynamic Stark shift, which has the 
same form as in the above case. This does not imply that there are no higher 
resonances at all, but rather that they are essentially suppressed. 
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3.2 Atom-classical field interaction 

We will proceed with the analysis of the interaction of atomic systems with 
classical fields without RWA. A generic interaction Hamiltonian has the form 

H = E UjXi + J2 9 3 (Xj + Xj) cos(n 3 -* + 0J), (36) 

33 3 

and admits some explicit resonances, u)j = ilj. The whole set of effective (im- 
plicit) interactions can be easily obtained in the same way as it was outlined in 
3.1. To be able to use the effective Hamiltonian in the form (|31|l we first rewrite 
the Hamiltonian (??) in the Floquet form by making use of the phase operators 
{Eqj, Ej}, which are generators of the Euclidean algebra: 

[E 0j , E k ] = -5 ]k E 3 , [E oj , Sj] = S jk E) , [e 3 , e(\ = 0. (37) 

Each group of operators (labeled with the same index j) acts in a Hubert space 
spanned on the eigenstates of the (Hermitian) operator Eoj : 

E 0j \n) j =n\n) h n = ...- 1,0,1,..., (38) 

so that in the basis (|38|l the operators Ej , Ej act as rising- lowering operators: 

E j \n) j = \n-l) j , E]\n)j = \n + l)j. 

The phase states the eigenstates of Ej, Ej\6)j = e" 11 ^ |z9)j, are not nor- 
malized. 

Now, let us consider the following time independent Hamiltonian 

H°» = E "fa + E + E 9i (*i + X]) {E j + E]), (39) 

j jj j 

It is easy to observe that the average value of the above Hamiltonian over the 
phase states in the rotating frame, 

|0(t)) = nj-e-^^l^, (40) 



OO 

— E 

N + 1 ^ 



\&)j = lim — 7 , 

v n— — oo 

where \n)j are the basis states 1|38[). coincides with the Hamiltonian (|36)1 due to 

(m\E 03 \m) = o, m)\Ej\m) = e-^ t+ ^\ 

The Hamiltonian (j39(l is a Floquet form of the initial time dependent Hamil- 
tonian (|36|l (from now on we put all the phases fij equal to zero). In the case 
of weak driven fields, g -C w, £1, the Hamiltonian (|39() can be represented in a 
form of expansion over the principal resonances which can be observed in this 
system according to (pTTJ) . 
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In the single channel case, the Dicke model in the classical field, we have the 
classical problem, first solved by Shirley |3] (for a single two level atom). Let us 
consider a collection of A two level atoms in a linearly polarized EM field. The 
corresponding Hamiltonian has the form 

H = wS z + g (S- + S+) cos fit, (41) 

where is the classical field frequency. 
The Floquet form of igTJ is 

H op = ljS z + QE + g (S- + S+) (E + E r ), (42) 

The following identifications 

X Q = S z , X± = S±, Y = E , Y+ = E\ F_ = E, 

so that M x o) = (1 + A/2)A/2 -X% + X Q , <j) 2 {Y a ) = E^E = EE^ = I, lead to 
Oh (S z , n) = 0, k > 2 which immediately leads to the effective Hamiltonian l|31|l 

H° p ff = loS z + nE + 2geE^ES z + g J iz^fL (E 2l+1 S+ + h.c) . (43) 

i=o L 

Averaging the above Hamiltonian over the phase states we obtain an effective 
time-dependent Hamiltonian 

H eff = cjS z + 2geS z + g £ t^t ^(2i + i)nt s+ + h ^ 

1=0 

which means that only odd resonances, f2 sa (21 + l)u>, appear in this system. 
This result is different from the atom-quantum field interaction (and surpris- 
ingly can not be obtained from the corresponding quantum Hamiltonian 135fl by 
substituting the field operators by c— numbers), because in this case the struc- 
tural function l|49|) for the Euclidean algebra (|37(l is a constant, that is, a zero-th 
degree polynomial function. 

4 Conclusions 

Quantum systems possessing as the integral of motion the operator correspond- 
ing to the total number of excitations (classically integrable systems) admit a 
simple classification of possible resonant transitions, which are separated into 
explicit (which appear in the original Hamiltonians) and implicit (effective) res- 
onances. The total number of such resonances is always finite and depends on 
the number of atoms and of atomic levels. By a simple algebraic manipulation 
of the free Hamiltonian and the integrals of motion, one can obtain some spe- 
cific conditions for the frequencies of interacting subsystems, which define the 
implicit resonances. Nevertheless, the form of the interaction terms associated 
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with each implicit resonance condition can be found only after obtaining the 
effective Hamiltonian (in a perturbative way). 

Quantum systems not preserving the total number of excitations (classi- 
cally non-integrable systems) do not admit a simple classification of admissible 
resonant transitions. The explicit form of the effective Hamiltonian is needed 
to specify the allowed resonances, and their number is always infinite. Such 
non-linear (on the generators of some Lie algebra) Hamiltonians can be repre- 
sented as a series of operators describing all the possible transitions, which might 
become resonant under specific relations between frequencies of interacting sub- 
systems. The structure of the effective Hamiltonian essentially depends on the 
algebraic structure of interacting subsystems (polynomials Px (X ) ,Py (%))• 
In particular, it is reflected in the types of resonances which are allowed for a 
given system. 

In the vicinity of each resonant transition all of the other transitions can 
be considered as a perturbation. Thus, in the course of evolution some spe- 
cific finite-dimensional subspaces in the Hilbert space of the whole system are 
approximately preserved, and for each of these invariant subspaces there is a 
corresponding integral of motion. 

Acknowledgement 1 This work is partially supported by the Grant ^5704 of 
ConsejoNacional de Ciencia y Tecnologia (CONACyT). 

5 Appendix 

The method of Lie-type transformation (small rotations) H5| provides a 
regular procedure for obtaining approximate Hamiltonians describing the effec- 
tive dynamics of nonlinear quantum systems. The idea of this method is based 
on the observation that several quantum optical Hamiltonians can be written in 
terms of polynomially deformed algebras sl pd (2,R) [Tfi ] . [17 ] . [T3 ] . |T3 ] . PH ] , 

H izlt =AX +g(X + + X.), (44) 

where the operators X± and Xo are generators of the deformed algebra and 
satisfy the following commutation relations (J5J, and 

[X+,X-]=P(X ), (45) 

where P(Xq) is a polynomial function of the diagonal operator Xq with coeffi- 
cients that may depend on some integrals of motion Nj . If P(Xq) is a linear func- 
tion of Xq, the usual sl(2, R) or su(2) algebras are restored. If for some physical 
reason (depending on the particular model under consideration) i) = j/A<l 
is a small parameter, the Hamiltonian (|44|l is almost diagonal in the basis of 
the eigenstates of X and can be approximately diagonalized by applying in a 
perturbative manner the following unitary transformation (a small nonlinear 
rotation) 

U = exp[r)(X + -X_)}. (46) 
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Applying the transformation l|4t)|l to the Hamiltonian 144fl according to the 
standard expansion 

00 k 

e A Be- A = Y J \^ k AB), (47) 

fc=i 

where adyt is the adjoint operator defined as &&a{B) = [A, B], we obtain 

oc , 

H cS = UH hn U ] = A Xo+gY ?? fc ■ -- ad£(X + + (48) 

where T = X + — X- and we have taken into account that, due to l|45|) . 
ad T (X ) = [T, X ] =-(X++X-). 

The effective Hamiltonian acquires the following form 

00 

H cS = AX + gvY'v* [Xlf k (X , v) + h.c] , 

fc=0 

where fk(Xo, rj) is a function of the diagonal operator Xq and can be represented 
as a series on r\: 

fk(X , V ) = i^±llv fc+1 0(^o) +0(rj), 

and 

0(X O ) = X+X- (49) 

is a structural function, V4>(Xq) = 4>(Xq) — 4>(Xq + 1) = P(Xq); the prime ( ') 
in the above sum means that the term with fc — is taken with the coefficient 
1/2. 

By keeping terms up to order 77 we get 

H cS = AX o + ry.gV0(X o ), (50) 

and in the first approximation the resulting effective Hamiltonian is diagonal on 
the basis of eigenstates of X - 

The higher-order contributions always have the form X+fk(Xo) + h.c. + 
g{Xo). This makes the procedure of removing the off-diagonal terms somehow 
trivial at each step, in the sense that it is always obvious which transformation 
should be applied. For example, to eliminate the terms of the form 

r, l [XXf k {X ) + f k {X )X% 

it suffices to apply the transformation 

cxpi^-Tk), (51) 

with Tfc = X+fk(Xo) — fk(Xo)X_, since the first commutator of A Xq with Tj. 
cancels the corresponding term in the Hamiltonian. 
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